Quantum Computation Based on Retarded 
and Advanced Propagation 

Giuseppe Castagnoli* 
Computation is currently seen as a forward propagator that evolves (re- 
tards) a completely defined initial vector into a corresponding final vector. 
. Initial and final vectors map the (logical) input and output of a reversible 

Boolean network respectively, whereas forward propagation maps a one-way 

> : 

O ■ propagation of logical implication, from input to output. 

z. 

Conversely, hard — NP-complete — problems are characterized by a two- 
way propagation of logical implication from input to output and vice versa, 

in ; 

given that both are partly defined from the beginning. Logical implication 

ON 



can be propagated forward and backward in a computation by constructing 
the gate array corresponding to the entire reversible Boolean network and 

o 

constraining output bits as well as input bits. The possibility of modeling 

^ : 

the physical process undergone by such a network by using a retarded and 

-i— > 

^ . advanced in time propagation scheme is investigated. 

3 . 

PACS numbers: 89.70.+C, 02.50.-r, 03.65.-w, 89.80.+h. 

> 

H 1 I. THE NOTION OF COMPUTATION REVERSIBILITY 

This work deals with the possibility of deepening the notion of computation reversibility 
in the quantum framework. This might provide a richer form of computation. 

The discovery that computation can be reversible in principle and the development of 
the thought model of a classical reversible computer, became cornerstones of research on the 
physical basis of computation [1, 2, 3, 4]. Once found to be reversible, computation could 
be moved into the quantum framework [5, 6, 7, 8, 9]. 
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The finding that the deeper the physical basis the richer is the form of computation, has 
been a remarkable result: 

• a quantum register undergoing a unitary evolution can perform in quantum super- 
position (or parallelism) a number of computations which grows exponentially with 
register size, 

• superposition elements interfere and the successive measurement brings about a col- 
lective property of those elements in the episystem, 

• this enables faster than classical computation in some special instances [7, 10, 11, 12, 
13, 14]. 

Although essential to bring computation into the quantum domain, the notion of com- 
putation reversibility entered the quantum framework unaltered. This, we argue, might 
constitute an undue restriction. 

As is well known, the notion of reversibility, thus of the absence of a privileged direction 
of causality, can undergo a quantum change. Models that postulate the non-redundant 
coexistence of a retarded and an advanced propagation in the quantum framework have 
been formulated [15, 16, 17, 18, 19], although mostly outside a computational context. 

FIG. 1. 

Let us examine the notion of computation under a broader perspective, keeping in mind 
that it can be mapped on either a one or a two-way propagation. Currently, quantum 
computation is seen as a propagator that evolves (retards) a completely defined initial vector 
into a final vector, which is therefore a function of the initial vector. One-way physical 
propagation maps the one-way propagation of logical implication from the inputs to the 
outputs of a reversible Boolean network (Fig. |Tja) . Such a network appears in the space-time 
diagram of a "sequential computation" process: network nodes are the register qubits before 
and after each transformation, Boolean gates and wires represent the transformations. 
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The existence of a privileged input-output direction (of one-way propagation) in sequen- 
tial computation sets some constraints on the underlying Boolean network. A Boolean 
constant must be preassigned to any input node, whereas none is allowed on output nodes. 
An output cannot be connected to an upstream input which participates in the logical de- 
termination of that same output (Fig. |ï]b; ignore dotted arrow). The fact that sequential 
computation is suspected to be unable to solve NP-complete problems in polynomial time 
[20] might be related to the presence of such constraints. 

Conversely, Boolean networks do appear in problem solving in a less constrained way 
An example is the problem of finding the satisfiability of a genèric Boolean network. Fig. [T]c 
exemplifies the Standard way of stating this problem for a reversible network. Out of hàbit, 
the two extremities of the network are still called input and output, although there is no 
longer a privileged input-output direction. Both the input and the output must be partly 
constrained: constant Boolean vàlues must be assigned to part of the input nodes, another 
part must be left free. The same applies to the output. The Boolean constants applied to 
both the input and the output propagate a conditional logical implication, respectively in the 
input-output and output-input directions (of course the logical implication, say, in the input- 
output direction is conditioned by the possible Boolean assignments of the unconstrained 
part of the input, and similarly for the other direction). The coexistence of two directions 
of conditional implication makes the problem of determining network satisfiability hard, 
namely NP-complete. 

The case of the Boolean network in Fig. (ignore the time axis) is similar. The 
output-input connection introduces another form of coexistence of a forward and backward 
conditional implication (along a loop), which makes the network satisfiability test equally 
hard. 

It is physically quite possible to propagate a conditional logical implication both forward 
and backward along a reversible Boolean network, by laying out the entire network spatially 
and by constraining part of the inputs as well as part of the outputs (Fig. |T]d). Thus 
network nodes r,s, ... become coexisting qubits |xi) r , \X2) S , (with Xi — 0,1) as time is 
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now perpendicular to the network layout. 

We will develop a thought model of this network by applying a sort of divide-and-rule 
strategy, meant to take advantage of a quantum symmetry related to particle statistics[|. 

Conceptually, gates, with their input and output qubits, are first thought of as being 
"divided" , all the constraints imposed on them through the wires are removed. Each inde- 
pendent (unwired) gate is implemented by a quantum transposition of simulated annealing — 
namely ground state computation. At equilibrium, each gate sits in a degenerate ground 
state. This can host any linear combination of those gate eigenstates — each a tensor product 
of the gate qubits eigenstates — which satisfy the gate Boolean relation. Disregarding wires, 
all gates would relax toward ground in parallel and independently of each other. This, at 
equilibrium, would produce a network state consisting of a superposition of tensor products 
which satisfy the Boolean relations imposed by the gates, while those imposed by the wires 
are ignored. 

Wires are restored by using a special symmetry induced by particle statistics. With no 

loss of generality, each wire is assumed to comprise a NOT function. As such, it will be called 

r s 

link. The truth table of link, say, r, s is thus 1 • Such a relation between the eigenvalues of 

1 

any two interlinked qubits r, s, will be imposed by a special symmetry A rs , related to particle 
statistics [defined by A TS \x) r |1 - x) s = \x) r I 1 ~ x) s , A rs \x) r \x) a = °> with X = 0, 1— see 
Section II. C]. Such symmetries, or their product for all links A = U r s A rs (the A rs pairwise 
commute), "rule" by restoring consistent wiring and network completeness. 

This kind of wiring is liable of two interpretations, one holding in the conventional one- 
way propagation model, the other in a two-way propagation model (by "propagation" , when 



So far quantum computation has used special features like superposition and multiparticle inter- 
ference, certainly not particle statistics. The possibility of exploiting particle statistics is investi- 
gated in this work. 
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not explicitly referred to conditional logical implication, we will mean the physical process). 
Both will be pursued throughout this work. 

In one-way propagation, one obtains a quantum annealing computation over a Boolean 
network. However, no advantage will be found over classical annealing (disregarding quan- 
tum tunnelling, not addressed in this work). 

In two-way propagation, it is as if A induced a watchdog effect on the overall network 
state \ip (t)) generated by the relaxation process of the independent gates. By the way, we 
should note that \tp (t)) is not necessarily a pure state, we will use the method of random 
phases. Symmetrization under A would operate as projection, canceling so to speak the 
amplitudes of those network eigenstates which would violate A, and reinforcing the others 
through renormalization. This would yield (roughly) to a population of gates relaxing in- 
dependently of each other, without ever incurring frustration, with an essential advantage 
over classical annealing^. 

In this scheme, any particle statistics symmetry S (for example, singlet or triplet state 
symmetry) must be seen as a projector S applying to the evolving state for all times t: 
Vi : S |t/> (í)) = \ifj (i)). In this view, S does not need to be an initial condition conserved as 
a constant of motion. Eventually, S commutes with the propagator of \ip (0)) ; but this can 
be so in a reversed perspective, where the propagator could be forged by projection. 

Considering a particle statistics symmetry S as the result of projection, is consistent with 
the usual interpretation where symmetry is seen as a conserved initial condition. However, 
the two interpretations are not necessarily equivalent: in the current computation context, 
the projection view appears to favour a two-way propagation model. We mean those inter- 
pretations of quantum mechanics which postulate that propagations are driven by both their 

2 To some extent, this is similar to present-day quantum error correction, where the onset of errors 
and decoherence is countered by repeated projection of a redundant computation state onto some 
symmetric subspace [21, 22, 23, 24]. 
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initial and final conditionst 15 ' 16 ' 17 ' 18 ' 19,33 l 

In current assumptions we are dealing with a set of independently relaxing gates, whose 
overall state, for all times t, is projected on the Hilbert subspace which satisfies symmetry 
A: propagation would thus be "driven" by relaxation and "forged" by projection. 

Let us consider first the usual one-way propagation model. Here propagations are mod- 
eled by calculus, which is of course one-way. Any infinitesimal relaxation should be submitted 
to projection, and so on repeatedly. However, this yields to a Zeno effect that "freezes" the 
propagation in its initial state: Ví : \ip (t)) = \ip (0)) (Section II. D). 

It might be argued that one-way calculus is not fit to model a propagation "forged" 
by projection. Intuitively, the nature of projection appears to be discontinuous; even if 
projection is applied at all times t, this does not necessarily mean that it can be used in a 
differential way As a matter of fact, the Zeno effect disappears when this projection-forged 
propagation is modeled by the two-way propagation scheme developed in [18]. 

Here the propagation of (i)) turns out to be driven by the final condition that pro- 
jection at all times t yields a state that simultaneously (under a system of mathematically 
simultaneous equations) approaches ground, according to the independent gate relaxation 
process, and satisfies A. 

Interestingly, projection on a predetermined outcome implies an elusive chronology vi- 
olation. At any time ti the propagation has to choose between two possible projection 
outcomes, one satisfying, the other not satisfying A. Satisfaction of A is the condition of 
that choice. The propagation meets this condition immediately after the choice: the choice 
should be before projection (at time íi), and the condition is met after projection (at time 
í 2 , reasonably t 2 > íi). Chronology violation is elusive since it can be confined within a time 
interval At = t 2 — íi as short as desired. It can be, so to speak, "safely" confined within a 
time interval At such as AtAE = fi, where the ordering of time and causality becomes fuzzy. 
However, the consequences of choices taken in Aí can naturally be "long range". For the 
sake of clarity, we should note that chronology violation would become completely evident 
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if the vector to be projected on a predetermined outcome were discretely different from it. 
A condition set in the future would determine a choice in the past. It becomes elusive when 
projection is performed for all í. 

Section II provides a thought model based on the above notions. It is not an imple- 
mentation model; some Hamiltonians are represented by Hermitean matrices. However, it 
appears to be plausible in principle. If implementable, it could speed up the solution of an 
NP-complete problem with respect to known classical computation. 

An analogy can be observed between this work and D. Deutsch's work on closed time-like 
lines [26] . Their relation to general relativity is irrelevant in this context — what is relevant is 
the presence of loops of conditional implication. In Fig.lb, forward and backward conditional 
implications coexist on the two halves of a loop which (in Deutsch's model) becomes a closed 
time-like line in a chronology violating region (consider now the dotted arrow). Sequential 
computation circulates along this line. This enables, without necessarily implying it, faster 
than classical computation (the problem was finding the fixed points of a hard to invert 
function). 

II. THOUGHT MODEL 

A. Classical Simulated Annealing 

With a view to introducing quantum annealing, we shall review a simple application of 
classical simulated annealing — also called ground state computation — to the satisfiability 
problem. Network nodes (Fig. ld) can be seen as cups, each containing a coin in either 
head or tail state (say Boolean or 1). Computation follows an if/then process: 

(1) does the coin configuration satisfy all the local relations imposed by gates and links? 

(2) If YES, stop: a solution has been found and the network is satisfiable. 

(3) If NOT, "shake" (randomly flip coins) and go to (1). 



7 



Network unsatisfiability is ascertained with any desired confidence level by looping for a 
sufficiently large number of times, without ever finding a solution. 

In order to speed up computation, each network element (N.E.) should be viewed as a 
mechanism whereby all states which satisfy the N.E. belong to a degenerate ground state, 
while the energy level of the others is discretely above ground (Fig. |I]d). This feature — 
during "shaking" — provides an energy gradient, orienting the evolution towards ground 
state. Dumping is of course needed. A heat bath of suitable temperature can provide for 
both dumping and shaking. 

Shaking is required because this kind of computation is subject to "frustration" . A coin 
configuration can satisfy most N.E. and still be far from any solution. The process, so 
to speak, should be undone. This boils down to the presence of relative energy minima — 
"shaking" is needed to avoid state trapping in such minima. Let N be the number of network 
nodes. The number of relative minima and relaxation time — expected time to reach the 
absolute minimum — is exponential in N, according to today's knowledge by any classical 
simulated annealing strategy. For a description of efficient strategies, see [27]. 

B. Quantum Annealing — Gates 

In quantum annealing, the Hamiltonian of the Boolean gate operating on qubits, say, 
s,t,u,v (Fig. |ï]d), has 16 eigenstates, namely all possible tensor products of the gate qubit 
eigenstates. Those tensor products which satisfy the gate belong to a degenerate ground 
state, while those which do not satisfy it are metastable states belonging to discretely higher, 
say by more than Ai?, energy levels. An external input or output with a preassigned Boolean 
value becomes a one-qubit gate, and that value will correspond to the ground state of the 
gate Hamiltonian (Fig. |T]d). 

Given a heat bath temperature T such that < kT << AE, in a semiclassical picture 
the probability p of finding (as the result of a measurement) the gate in ground state by the 
time t follows asymptotically the law 
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p = 1 - e-*, (1) 

where a > 0. Of course this exponential law holds when relaxation has reached a constant 
rate. Model Hamiltonians of elementary Boolean gates with coexisting inputs and outputs 
can be found in [28]. Of course such gates differ from the gates used in the sequential 
approach, where inputs and outputs are successive states of the same register [6, 7, 29, 30, 
31, 32]. 

Links could be implemented accordingly. Then, link-gate interaction Hamiltonians 
should be added. The number of local minima and computation time would be likely to re- 
main exponential in N (the benefit deriving from quantum tunnelling through local minima 
is not considered in this work). Links will be dealt with in a different way 

C. Quantum Annealing and Particle Statistics — Links 

The correlation between simultaneous qubit eigenvalues imposed by links will rely on 
a special particle statistics symmetry — seen as a projection operator (Section I). The 
(thought) physical situation giving origin to such symmetry (one for each link) is still ob- 
tained by means of a relaxation process. 

Without any restriction, a network topology is adopted where nodes come in interlinked 
pairs (Fig. [!]). The "divide" rule is still applied. By "turning off" gate Hamiltonians, links 
are, so to speak, "ungated" . This leaves us with a set of independent links, which comprise 
all the nodes of the network. 

For the sake of visualization, each link node r, s (Fig. |]d) can be seen as a site of a spatial 
lattice. These sites are labeled by A = r, s. Each site hosts exactly one spin | particle, with 
its a z spin component either down (x = 0) or U P (x = !)■ Therefore qubit \x)\ reads: the 
o z component of the particle hosted by site À is x- However, we should like to keep this 
model abstract: À and x can be any two binary compatible attributes (dynamical variables) 
of a fermion. 
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In a particle label representation, qubit \x)\ is written \x) 1 |À) 1 ; for example |0) r becomes 
|0) 1 |r) 1 . By this we mean that the x attribute of particle 1 has eigenvalue and that its 
À attribute has eigenvalue r. Or, since we are dealing with two identical partides, particle 
1 is interchangeable with particle 2; in this case |0) r becomes |0) 2 |r) 2 . Therefore, qubit 
representation will be the outcome of a finer identical particle representation. 

The lattice L of network nodes (or qubits), with exactly one fermion per node (one 
independent \ attribute per qubit), is created as follows. Gate Hamiltonians are still "turned 
off" . Link Hamiltonians, thus completely independent of each other, are "turned on" . Then 
links are allowed to relax on their ground states. 

The Hamiltonian H x 2 of link r, s is defined as follows. It operates on the spatial coor- 
dinates of fermions 1 and 2 — not on their spin. Therefore H x 2 operates on the Hilbert 
space: 

f 1 1 

H x = span | \r) 1 \r) 2 , \s) 1 \s) 2 , -j= (\r) 1 \s) 2 - \s) 1 |r) 2 ) , -j= (\r) 1 \s) 2 + \s) 1 \r) 2 ) 

where |r) 1 |r) 2 reads particle 1 in site r, particle 2 also in site r, etc. 
The energy eigenvalues <-> eigenstates of the Hamiltonian 

^ tp\ n n n \ 



H 12 



E x 








E x 
















E x <-> | 




E x <-> | 



, with E x discretely above 0, are 



r) 1 \r). 



^ ^ (k> a |«> 2 - |*> x |t-> 2 ) 
<-> ^(\r) 1 \s) 2 + |s) 1 |r) 2 ). 

When link r, s has relaxed onto the ground state of H x 2 , each link node or qubit is associated 
with exactly one particle. This creates the "lattice" of qubits r and s. For the sake of 



10 



simplicity, the two lowest eigenvalues are both 0; they could as well be differentf^ provided 
they are well below E x . The link ground state in H x \ = H x ® H\ has thus four degenerate 
eigenstates (the overall state must be antisymmetrical): 

0^^|0>i|0> a (|r>i|a) a -|a)i|r> 2 ) 

0^^|l>i|l) a (|r>i|a> a -|a)i|r> a ) 

0^3 (|0>! |l> a + |1> X |0> 3 ) (|r> x |*> a - |t-> 3 ) 

- ± (|0> Jl> 2 H1>i |0> 2 ) (|r> Js> 2 + |s> Jr>^ 

To avoid triplet states of parallel spins, we introduce the Hamiltonian H x 2 operating on 
the Hilbert space 

H x = spanjlO)! |0) 2 , |1) 2 , -L (|0> 1 11) 2 + (l), |0) 2 ) , -L (|0> 1 11) 2 - (1), |0) 2 )| , 

where |0) 1 |0) 2 reads: the a z component of particle 1 is down, the same for particle 2, etc. 

( E* 0^ 



n \2 



E x 





with E x discretely above 0. 



H x 2 commutes with H x 2 and lifts the energy level of the triplet states of parallel spins. The 
link ground state remains two-fold degenerate and hosts the spin singlet state \if))_ and the 
antiparallel spin triplet state |^>) , : 



m T = 2 d°)i |l>a T I 1 ): I°) 2 ) GOi l s ) 2 ± l s >i l r ) 2 ) , 
upper (lower) signs go together. With simple manipulations, (2) can be written: 



(2) 



V2 
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(3) 



where 



3 which would be necessary if we regarded H x 2 as acting on spatial variables 
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W) = ^ (10): \r) x |1> 2 \s) 2 - II)! \s) l |0) 2 |r) 2 ) = |0> r , (4) 

V>") = ^= (|1>! |r>! |0) 2 |s> 2 - 10), 1^ |1) 2 |r) 2 ) = e « |l) r |0), . (5) 

It is readily checked that both terms appearing in in the particle representat ion 

map onto |0) r |l) s (|l) r |0)J in qubit representation; 5 is an undefined phase. Thus 

l^ = ^(|0)Jl) s T|l) r |0)J, (6) 

where we have set S = for \ip) to behave for a rotation of the x reference in eq. (6) as in 
eq. (2). 

Since \ip)_ and are degenerate, we can say as well that the two orthogonal ground 
eigenstates of link r, s are |0) r |l) s and |l) r |0) s . Therefore the degenerate link ground state 
has the form: 

ie) = «|0) r |l) s + /3|l) r |0) s , (7) 

with \a\ 2 + \(3\ 2 = l. 

Let A12 = 1 — P12 followed by renormalization be the antisymmetric operator in par- 
ticle representation. By definition, the application of A í2 to a genèric state filters out 
and normalizes its antisymmetrical component. We can see that A 12 \ip') = \ip') and 
A l2 \ïp") = |V">- In equivalent terms: A 12 |0) r |l) s = |0) r |l) s and Ai 2 |l) r |0> s = |l) r |0) s . 
Also: A l2 |0) r |0) s = |0) r |0) s and A 12 \l) r |l) s = |l) r |l) s , if we are dealing with the excited 
states of energy eigenvalue E x . However, such states are forbidden in the link ground state. 
They should not be confused with the "bad" tensor products: 

|0)J0) S =-L|0) 1 |0) 2 (|r) 1 | S ) 2 + | S ) 1 |r) 2 ), 

|1)J1) S =^|l) 1 |l) 2 (|r) 1 | S ) 2 + | S ) 1 |r) 2 ), 

which are symmetrical and violate A 12 . We should note that qubit notation, and qubit 
density matrices (|0) r (0| r , etc), on which the relaxation process operates, are the same for 
both kinds of states. 
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All this becomes, in H rs = span {|0) r |0) s , |0) r |l) s , |l) r |0) s , |l) r |l) s }, the projector (or 
"symmetry" ) A rs defined by the following equations: 

A-.lx> r |i-x>. = lx> r |i-x>., As\x) r \x) s = o, (8) 

with x — 0,1. A rs symmetrization projects a genèric state of H rs onto the symmetrical 
subspace Ha ts spanned by the legitimate link states, given by eq. (7). It is understood that 
the application of A rs is followed by renormalization. Let (L) be the expected value of the 
link energy. Clearly, A rs projection is equivalent to the two mathematically simultaneous 
conditions: A í2 projection and (L) = 0. 

D. A divergence between one and two-way propagation models 

The link state, in the reference where H* 2 is diagonal, is given by eq. (7). Let us con- 
sider continuous transformations of state (7) which leave A rs symmetry unaltered. Without 
significant restriction, state (7) is rewritten: 

|£>=costf|0>Jl) s + smtf|l>jO> s . (9) 

There is of course a unitary operator R(<p), acting on the Hilbert space H rs = 
span{\0) r \0) s , |0) r |l) s , |l) r |0) a , |l) r which changes |Ç) of eq. (9) into: 

10 = cos (0 + <p) |0) r |1), + sin (0 + <p) \l) r |0) s , (10) 

without ever (for all ip) violating A rs symmetry; here ip = ut. Since •& and (p will be used to 
model the result of gate relaxation, we can assume <p ^ and ^ -& + <p ^ |. 

It could be seen that R (ip) operates on the overall link state in an irreducible way 
However, for reasons that will become clear, we are interested in a factorizable operator: 

R rs { í p) = R r (y) R s i 1 ^), where R r (ip) [R s (ip)] operates on the state of qubit r (s), indepen- 

( ■ \ 

cos ip —smip 

dently of the state of qubit s (r). Let R t (<p) = , with t — r,s. 

sin ip cos ip j 

We should note that R rs (p) does not commute with A rs . Its application to state (9) 
introduces the terms |0) r |0) s and |l) r |l) s . There are two possible cases: 
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(i) these terms are the excited states which satisfy A12, then the link goes into a super- 
position of two energy eigenstates, ground and E x , while its expected energy (L) is 
raised, and 

(ii) these terms are the symmetrical states of paral·lel spin, which violate A12] such states 
would be, so to speak, immediately projected off. 

In this latter case R rs (tp) operates adiabatically on the link, leaving it in its ground eigen- 
state, where A rs is satisfiedQ. Further below, it will become clear that this is a possible 
solution. 

We shall formalize this, first in the conventional one-way propagation model. A series 
of infinitesimal operations R rs (Sp), each followed by A rs symmetrization should be applied 
until f óp = p. However this operation suffers from a Zeno effect which keeps state (9) 
unaltered: indeed lim n _,.oo [A rs R rs |£) = as it can be checked. In fact, in first order, 

R rs (tp) drives the link state out of the symmetrical subspace Ha ts = span {|0) r |l) s , |l) r |0) s } 
orthogonally to it. A rs projection on Ha tb brings the link state back to the starting point. 
Therefore R rs (<p) cannot operate adiabatically in a one-way propagation perspective, or in 
calculus which is naturally one-way. It would inevitably lift (L). 

However, calculus may not be applicable to a situation which involves a continuous 
projection due to A rs symmetrization. As a matter of fact, during an elementary rotation 
Aip = 2 , symmetry A rs is continuously violated (if the link remains in its ground state, the 
fundamental symmetry A12 is violated). This error may not vanish in the limit Aip — > 0. 

A reasonable alternative is to think that rotation and symmetrization should be applied 
in a mathematically simultaneous way, so that A rs is never violated. This removes the Zeno 
effect and brings in a two- way propagation (Section II. F). We consider a free normalized 
vector |£ ^ of the Hilbert space H rs . By this we mean that its amplitudes on H rs basis 
vectors are free variables (up to normalization) independent of each other. This vector 

4 For a study of adiabatic operation on energy eigenstates, see [33]. 
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must satisfy, for all ip or t (<p = ut) , the following system of mathematically simultaneous 
equations: 

(i) As\0= IO> 

(ii) Tr s |f) = cos 2 (0 + <p) |0) p (0| r + sin 2 (0 + |l) r (l| r , 

(iii) Tr r |f) <Ç"| = sin 2 (tf + |0} s (0|, + cos 2 (0 + <p) \l) s (l\ s . 

Condition (i) ensures that A rs is never violated and is equivalent to A 12 £') = and 
(L) = 0. Conditions (ii) and (iii) mean that the density matrix of each qubit has been 
rotated by <p. For all (p, the above system of mathematically simulteneous equations yields 
the solution: 

iS , 



cos (t? + tp) |0) r |l) s + e là sin (0 + ^) |l) r |0) s , 

as it is readily checked. The further condition that || (£| Ç") || be maximum (see Section II. F), 

sets 5 = 0, thus yielding to state (10). 

As we will see (Section II. F), this is exactly an application of a two-way propagation 

modelQto the projection process induced by A rs (by A 12 under an adiabatic operation). To 

sum up, in this scheme, one can operate separately and adiabatically on the link qubits; the 

link state undergoes a unitary evolution while remaining in its degenerate ground state, thus 

always satisfying A rs . This unitary evolution is dríven by conditions (ii), (iii) and forged 

by condition (i), namely continuous A rs projection. The fact that the resulting evolution is 

unitary thus reversible (it does not dissipate free energy), is consistent with the assumption 

of adiabatic operation. 

We will show another way of obtaining the same result. Let qubit r be connected to the 

E r 0^ 



one-qubit Hamiltonian 



V 00 / 



operating on span {|0) r , |l) r } and qubit s be connected 



5 We will use the model developed in [18, 19], whose outline will be given in Section II. F. 
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to the one-qubit Hamiltonian 



yOE. ; 



operating on span {\0) s , |l) s } . Disregarding A rs , 



the two qubits would relax independently of each other onto |l) r |0) s (a heat bath should 
provi de for both "shaking" and "dumping" ) . 

The expected energy of the two gates in state (9) (the initial condition) is E 1 = 
(E r + E s ) cos 2 1?. Now we set the final condition that the state of this network simulta- 
neously satisfies A rs and has a lower expected energy, say E 2 = (E r + E s ) cos 2 ($ + </?). We 
assume that the energy difference has gone to the heat bath. This yields to state (10) again, 
as it can be checked. Now the unitary evolution from (9) to (10) is driven by the independent 
gate relaxation process, and of course it is forged by condition (i). 

One can think that gate relaxation may not always operate adiabatically on the link. The 
excited (antisymmetrical) components |0) r |0) s and |l) r |l) s would appear in a superposition 
with the link ground state. However, assuming E r ,E s <C E x , these components would 
have a small amplitude and should consequently be short lived, due to the continuous gate 
relaxation. The above two-way evolution would still apply to the Everett worlds where the 
link is in its symmetric ground eigenstate. 

We should further note that, rotating either independent qubit is equivalent to rotating 
both; the two rotations are redundant with respect to one another. This can be checked by 
removing either condition (ii) or (iii); just one is suflicient to yield state (10). 

We shall now go back to the one-way model. It can be checked that there is no way of 
operating separately and adiabatically on the two qubits to perform a unitary transformation 
of the link state inside the Hilbert space H Ars . Gate relaxation will raise the link expected 
energy (L) by some amount, non-infinitesimal this time. Therefore the link state can undergo 
a rotation. Then the link would relax on ground restoring A rs on the rotated state. This 
would go on repeatedly. However, this form of quantum annealing works exactly like classical 
annealing. 
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E. Turning Gate Dynamics on 



To sum up, having turned link Hamiltonians on and allowed links to relax onto ground 
(independently of each other), the entire lattice of qubits L is obtained, where qubits eigen- 
values are submitted to the symmetry (affecting any pair of qubits connected by a link) 
A = ïl r ^ s A rs . Note that all A rs operate on disjoint Hilbert spaces, concerning the states of 
disjoint pairs of nodes, and are pairwise commuting. 

After creating L and A, gate Hamiltonians are "turned on", then gates are allowed to 
relax. Here we consider the network at equilibrium, when all gates and links are in ground 
state. How this state is reached will be the subject of Section II. F. 

Let P g be the projector on the ground state of gate # g, and = J2h a h \^h) the genèric 
network state, where \^h) is a tensor product of all qubit eigenstates. At equilibrium, when 
all gates are in ground state, for all P g each \^> h ) should satisfy P g \^ h ) = \^> h ) . This yields 
to | x I / / l )s which satisfy all gates. Each \*$> h ) should also satisfy Aj^) = \^h)- This yields 
to|^^)s which satisfy all links. All P g operate on disjoint sets of qubits, thus on disjoint 
Hilbert spaces, and are pairwise commuting. Let P = Ií g P g , where g runs over all gate 
labels. The foregoing simultaneous equations are summarized by: 

ÍP|tt) = |tt), (11) 

Any \^h) satisfying eq. (11) satisfies all gates and links and is thus a solution (for the case 
with no solution, see Section II. A). Substituting AP with PA yields to the same conclusion. 
Even if A and P did not commute, the network ground states are eigenstates of both. 
However, the propagator of the actual relaxation process is shaped by A and commutes 
with it. 

This approach would suffer no interaction with the environment provided L and A remain 
in place: entanglement with a reservoir initially in state \R) changes each \^h) \R) into 
\^h) \Rh), where |i4) is now correlated with \^h), without altering the selection performed 
by eq. (11). 
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F. Relaxation Time 



We shall now discuss relaxation time. Let us consider a simple "network" of just one link 



r, s and one qubit gate, whose Hamiltonian operates on span {\0) s , |l) s } and has the form 
E s o\ 

. Let the initial state of this network be the symmetrical state (ti)) = \0) r \ l) s . 

\° °), 

For the sake of simplicity, only qubit r is assumed to be in direct interaction with a heat 
bath. At equilibrium, the (symmetrical) ground state should instead be |l) r \0} s . 

We shall apply the one-way propagation model first. We will limit ourselves to examining 
the network initial evolution from |0) r |l) s . This is readily done and provides useful insight. 
Since the initial state is factorized, relaxation of qubit s and heat bath perturbation on qubit 
r are (initially and in first approximation) independent of each other. 

The initial independent relaxation of gate s can be modeled by the two-step evolution: 

I* (*i)>. = 1 1>. ^ I* = c °s ^ |1>. + Ay? |0) s , (12) 

where t 2 > íi, and 5 is an undefined phase due to entanglement with the outside world. 
We are using the method of random phases: the average over 5 of the density matrix 
corresponding to the state on the right of the arrow, yields cos 2 Aip \l) s (l| s +sin 2 Ay? |0) s (0| a , 
as expected. 

The initial independent heat bath perturbation on qubit r can be modeled by the two- 
step evolution: 

IV> (íi)) r = |0) r -> \^ (t 2 )) r = cos Atf |0) r + e 5 ' sin Atf |l) r , (13) 

where Aí9 is a random phase. 

We can approximate the network initial evolution by symmetrizing the tensor product 
of the two qubits: 

liK*i)> =|o> r |i> - - 

|^(í 2 )) =cosA^cosAi9|0) r |l) s + e l5 "sinAv9sinA^|l) r |0) s , (14) 

18 



up to normalization; here 5" = S + 5' . 

This is an interesting result. If A$ = 0, the evolution remains frozen in its initial state 
|0) r |l) s (see eq. 14). The relaxation of gate s is conditioned by the action of the heat bath 
on qubit r. This means that, so to speak, by a lucky coincidence the network can relax, 
otherwise it cannot. Likely, this would mean frustration in a more complex network. In this 
approximation, one can see no reason why quantum annealing should perform better than 
classical annealing. 

Before applying the two-way propagation model (TW) developed in [18, 19], we shall 
outline it. This model was meant to justify the projection undergone by a quantum state 
during various measurement operations. 

Both forward and backward propagations are referred to the same direction of time in 
the episystem, on the basis of the following argument. Time can be seen as a parameter 
in the endosystem. It is the direction of causality in the endosystem — either forward or 
backward in time — that makes the difference between the two propagations. Therefore, 
both propagations are represented by kets evolving forward in time: one associated with 
causality forward in time, the other with causality backward in time. Because of the absence 
of a privileged direction of causality, the two propagations should be indistinguishable. 

Indistinguishability requires that (i) each propagation undergoes the same transforma- 
tions of the conventional one-way propagation, up to an overall phase, which is irrelevant 
to each propagation, and (ii) the relative phase of the two propagations is either or n. 

As a consequence of indistinguishability, the two propagations should proceed in coherent 
superposition. This superposition is the complete description of the propagation. Both 
initial and final conditions, comprising normalization, should be applied to it, whereby the 
initial (or indifferently) final condition of each propagation is just a free un-normalized state 
(completely undefined before the application of all conditions) of the system Hilbert space. 

State projection (due to either a measurement or a watchdog effect) implies a n shift 
in the relative phase of the two propagations. Consequently, their superposition changes 
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from constructive to destructive interference or vi ce- versa. Clearly, this can justify the 
discontinuous evolution of the quantum state throughout a projection. 

The foregoing model yields to the following simple rules to be applied to the superposition 
(the complete description): 

(i) the state before reduction \ip (ti)) coincides with the usual state; 

(ii) the state after reduction \ip (í 2 )) is a free normalized state of the system Hilbert space, 
subject to the final conditions and to the supplementary condition that 

(iii) ||(^(ti)| (t 2 ))\\ is maximum. 

A preparation, or a propagation thereof, can be the state before projection. The final 
conditions can be set from the future, for example they can be the result of partial or 
complete measurement, unforeseeable from the past. 

As it should be clear now, in Section II. D we have applied the two-way propagation 
model schematized above. The result was that A rs works like interference (destructive due 
to projection, constructive due to renormalization) on the relaxation of independent gates. 

We shall now apply the TW model to the simple network introduced at the beginning 
of this Section. 

The independent relaxation of gate s is modeled by the two-step evolution (12). The 
gate energy in \ip (t 2 )) s is E s cos 2 A(p (from now on, by energy in a state we will mean the 
expected value of the energy operator in that state). 



represent the energy of the qubit-heat bath interaction. Its expected value in \ip (t 2 )) r is 
E r cos 2 Ad. A$ can be assumed to be a random variable with uniform distribution in 
[0,27r], consequently cos2Aï? is zero on the average. We further assume E r <C E s . 

Since the gate and the heat bath operate independently (given that A rs works as inter- 
ference), the network energy at time t 2 is E = E s cos 2 A<p + E r cos 2A$. 



( 




The heat bath perturbation on qubit r is modeled by (13). Let the operator 



V 



r 



20 



The network state at time t 2 is a free normalized state of the network Hilbert space, 
submitted to the final conditions: 

(i) |\P(í 2 )) satisfies A rs ; 

(ii) the energy in \ip (í 2 )) is E. 
This yields the two-step evolution 

IV (*i)> = |0> r |l) s - IV fe)) = cos A^' |0) r + e lS " sin Ay.' |l> r |0) s , (15) 

where Atp' is still to be determined. The former condition (iii) is irrelevant here. 

The energy in \ip (í 2 )) is -E s cos 2 Ay/ + E r cos 2A<p'. Therefore Ay/ must be a solution 
of the equation 

E s cos 2 Aip' + E r cos 2 A(p = E s cos 2 Ay> + E r cos 2 Ai?. 

It can be readily checked that there is a unique solution in the interval of interest. If 
E r <C E 8 , A(p' — i.e. network relaxation — closely follows the independent gate relaxation 
Aíp. Moreover, Aip = Aip' on the average. 

This is quite different from the one-way model. Since the network energy is the sum of 
the energies of the independent network elements, frustration is never met. For example, 
the absence of a heat bath (A-d or E r = 0) does not freeze relaxation, unlike in the one-way 
model. 

This can be generalized to any Boolean network. Condition (iii) , moving in time together 
with continuous projection, allocates the overall energy among the network degrees of free- 
dom. Still, the overall energy decreases as if all network elements relaxed independently of 
each other and frustration is never met. Although a dynamic model is still lacking, this 
qualitative difference should yield an essential advantage over classical simulated annealing. 

Let us note that evolution (15) can represent an identical rotation of the two independent 
qubits r and s, obtained by rotating just qubit s. We must put E r = (thus obtaining 
Aíp' = Aip) and 5" = 0, which ensures rotation additivity; this is legitimate since we are 
now dealing with a pure state rotation. We find again the result of Section II. D. 
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III. DISCUSSION 



We have developed a form of quantum annealing computation where network wiring is 
implemented by a quantum symmetry A related to particle statistics. How this quantum 
wiring operates can be interpreted in two different ways. 

In a one-way, conventional, propagation model, this wiring does not exhibit any special 
feature with respect to classical wiring. Disregarding tunnelling, this quantum model does 
not yield any improvement over classical annealing. 

In an unconventional — two-way — propagation model, particle statistics works as a 
projection. A would continuously project the state of the independently relaxing gates 
on a "symmetric" subspace where network wiring is satisfied. This is naturally a blunt 
conjecture — this work should be considered an exploration of the possible role played by 
particle statistics (viewed as projection) in quantum computation. 

We should remark that something similar appears in other situations involving particle 
statistics. For example, given two non-interacting identical spin | fermions whose spatial 
wave functions overlap, there is a correlation between their mutual distance and the character 
of their spin state, either singlet or triplet. This correlation, due to antisymmetrization, 
creates a sort of "wiring" between two otherwise independent parts. 

Although an implementation model is still lacking, we see no reason to think that this 
is out of reach. In the example where À labels the sites of a spatial lattice and x labels the 
a z component of the spin \ particle hosted by each site, the link behaviour is likely to be 
related to exchange interaction. 

Were this model implementable — for example in the case of the single link and a one 
qubit Hamiltonian (Section II. F) — relaxation speed could be tested. In principle, this 
should decide which propagation model — either one or two-way — applies. 

The hypothesis of a computation speed-up based on chronology violation links this work 
with D. Deutsch's paper on computation along closed time-like lines. From that paper 
we quote: "it is curious that the analysis of a physical situation which might well not 
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occur should yield so many insights in quantum theory. But that is the nature of thought 
experiments. Perhaps we should also bear in mind that a frequently observed effect of time 
is to convert thought experiments into real ones" . Placing chronology violation within the 
time-energy uncertainty relation and the projection of an evolving state induced by particle 
statistics, might be a step towards feasibility. 

This research has been supported by Elsag Bailey, a Finmeccanica company. Thanks are 
due to A. Ekert, D. Finkelstein, L. Levitin, S. Lloyd, D. Monti, M. Rasetti and T. Toffoli 
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